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^ Identifying the symmetry of the superconducting order parameter in the recently- 

O discovered ferro-oxypnictide family of superconductors, RFeAsOi_a^F^, where i? is a rare 

I earth, is a high priority. Many of the proposed order parameters have internal tt phase shifts, 

Qh like the d-wave order found in the cuprates, which would result in direction-dependent phase 

c/^ shifts in tunnelling. In dense polycrystalline samples, these phase shifts in turn would result 

^ in spontaneous orbital currents and magnetization in the superconducting state. We perform 

H scanning SQUID microscopy on a dense polycrystalline sample of NdFeAsO0.94F0.06 with 

F^ Tc = 48 K and find no such spontaneous currents, ruling out many of the proposed order 

G parameters. 

O 

. ^. KEYWORDS: ferro-oxypnictides, pnictides, scanning SQUID microscopy, order parameter, Nd- 

> 

^It The recently-discovered ferro-oxypnictide family of superconductors includes materials 

^^ with transition temperatures above 50 K-^ and shows evidence that competing magnetism 
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plays a key role in the superconductivity.^"^ Determining the superconducting order parameter 
(OP) is key to understanding the interactions that induce superconductivity, but the OP of the 
ferro-oxypnictide family of superconductors remains uncertain. We report a phase-sensitive 
k>i test of the symmetry of the OP using scanning magnetic microscopy of dense polycrystalline 

samples. 

Grain boundaries form naturally occuring Josephson junctions that can carry supercur- 
rents. It is now well-known that the OP of the cuprate superconductors contains tt phase shifts 
associated with the d-wave symmetry, and that a tt phase shift can result upon going around 
a closed path in a polycrystalline sample; whether there is a tt shift depends on the relative 
lattice and interface orientations of the grains along the loop,^ for example as diagrammed in 
Fig. 1. TT-loops result in orbital frustration and spontaneous currents, as demonstrated by ob- 
servation of half-integer flux quanta in tricrystal cuprate samples.^ In polycrystalline cuprate 
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Fig. 1. Spontaneous magnetic flux generated by various configurations of grains, (a) Three coupled 
grains of a d-wave superconductor with no orbital frustration. The orientation of the OP in each 
grain is indicated; shading indicates the sign of the OP and dashed lines directions of strong inter- 
grain tunnelling, (b) Coupled grains with orbital frustration. If Xj is smaller than the grain size 
a half- flux-quantum (<l>o/2) vortex is generated, which can be positive (light) or negative (dark), 
(c) A polycrystalline sample with several <l>o/2 vortices, which will tend to couple antiferromag- 
netically. 

samples there is a finite density of tt loops, which, in well-connected samples (i.e. with the 
intergrain Josephson penetration depth Xj comparable to or less than the grain size), results 
in complex patterns of magnetization.^ 

Proposals for the ferro-oxypnictide OP include extended-5 order, with a tt phase shift 
between the hole and electron Fermi sheets, ^"-'^^ dx2_y2,^'^^~^^ dxy^^ p^^^^^^ s + d^^'^^ and 
s + id.^^ Phase-sensitive tests of the OP are important: because of the multiple Fermi sheets, 
and because time-reversal symmetry-breaking (TRSB) OPs remain a possibility, establishing 
the presence or absence of nodes does not definitively settle the OP symmetry. 

Which OPs would result in orbital frustration in a polycrystalline sample? Any pure d 
order will result in a tt phase shift between a and b axis tunnelling, leading to frustration. 
In principle a TRSB component could reduce the degree of frustration: drj.2_y2 + idxy order 
on a radially symmetric band, and in a 2-D sample (the c-axes of the grains aligned), does 
not give frustration. However, with the electron pockets, the ferro-oxypnictide Fermi surface 
is not radially symmetric, and any such reduction would be minimal, p order would result in 
frustration with or without a TRSB component, vr-shifted s order in principle could result in 
a TT shift between a and c axis tunnelling, if different, yr-shifted, sections of the Fermi surface 
dominate a and c axis tunnelling. Whether this is likely requires calculation. At present we 
must assume that an absence of spontaneous moments does not rule out vr-shifted s order. 

We have performed scanning SQUID imaging of a polycrystalline sample of nominal com- 
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position NdFeAsO0.94F0.06 grown by a high-pressure synthesis method. ^^ The superconduct- 
ing transition onsets at 51 K, with a midpoint at 48 K and a 10%-90% width of 2.7 K. The 
grains are weh-coupled: magneto-optical imaging and remnant magnetization measurements 
on a sample from the same batch indicate a bulk critical current of ^ 2000 A/cm^ at 5 K.^^ 

Our SQUID is a niobium-based scanning susceptometer design. ^^ Fig. 2(a) contains an 
image of the front end of the SQUID; magnetic flux is coupled into the 4.6 /xm diameter 
pick-up coil (the inner coil), the leads to which are shielded. In this SQUID a signal of 1 
$0 = hc/2e — 2.07 x 10~^^ Wb corresponds to a mean Bz in the pick-up coil of ^ 0.125 mT. 
The larger loop around the pick-up coil is a field coil; a measure of the local susceptibility can 
be obtained by applying a local field with this coil and measuring the response in the pick-up 
coil. 

We polished the sample to a shiny surface using AI2O3 polishing paper without any lu- 
bricant. In order to allow comparison with the vacuum, we scanned an area going over the 
edge of the sample. Our main results, scans of the sample cooled in different fields, are shown 
in Fig. 2. At fields below ^ 2 jiT (panels (b)-(e)) individual vortices are clearly resolved. 
They appear in different places on cooling in different fields, indicating that they are not 
frustration-induced spontaneous moments. For confirmation that they are regular vortices, 
they can be integrated: after subtracting a planar background, the three vortices indicated 
in Fig. 2(c) integrate to, from top to bottom, 1.05, 1.00 and -1.06 $0 (with a 5% systematic 
uncertainty due to uncertainty in the effective pick-up coil area). 

The 4.4 K scans reveal other features: several surface dipoles, clusters of vortices hinting 
at lumps of a magnetic impurity phase, possibly beneath the surface, and a widespread mot- 
tled background. In the area indicated in Fig. 2(b), the root-mean-square amplitude of this 
background signal, after plane subtraction, is 1.7 m$o- By lifting the SQUID slightly above 
the sample the sample can be heated while maintaining the SQUID below its Tc; a T = 55 K 
scan (Fig. 2(f)) confirms the presence of magnetic impurity phase. The surface dipoles also 
persist at 55 K while the mottled background disappears at Tc- 

Fig. 2(g), a scan at ^ 100 /xT, makes clear the granular nature of the sample: vortices 
cluster strongly in areas of weaker superconductivity. 

Fig. 2(h) shows sections of an isolated vortex and surface dipole from Fig. 2(c). The surface 
dipole provides a measure of the achieved imaging resolution: the peaks are separated by 5 
jivn. Assuming a 4.6 /xm-diameter pick-up coil and a point-like dipole this indicates a scan 
height of the pick-up loop above the sample surface of ^ 3 /xm. 2-D fits to the smallest dipoles 
in Fig. 2(f) also indicate a scan height of ^ 3 /im. 

Fig. 2(i) is a histogram of the full- width half-maxima of the vortices in Fig. 2(c). The 
narrowest are ^ 8 /xm. With a 4.6 /xm SQUID at a scan height of 3 /xm a vortex in a sample 
with zero penetration depth would appear with a FWHM of 6.2 /xm; i.e. the 8 /im FWHMs 
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Fig. 2. (a) Optical image of the front end of the SQUID; the field coil is highlighted in orange, the 
pick-up coil in green, and the shield of the pick-up coil leads in blue, (b-g) Images of the magnetic 
fiux coupled into the pick-up coil at sample temperature and applied field (b) 4.4 K, -3.3 /iT; (c) 
4.4, -2.7; (d) 4.4, -2.2; (e) 4.4, 0; (f) 55, 0; and (g) 4.4 K, +98 /iT. An ambient z-axis field of 
^ 3.3 /iT adds to these fields. The area indicated in (b) is used to analyze the background and 
the areas in (c) are integration areas (see text), (h) Cross sections along the green lines in (c): one 
section of a surface dipole and two of a vortex, (i) histogram of observed FWHMs of the isolated 
vortices in (c); two FWHMs per vortex, along the narrow and long axes. Hatched area indicates 
the approximate resolution limit. 

are strongly resolution-limited. However most of the vortices have observed FWHMs in the 
range 10-16 /xm, and most have visibly irregular shapes, suggesting that the actual vortices 
in the sample are spread out, with widths in the range of microns. 

A susceptibility scan (Fig. 3) shows which areas are superconducting: over these areas the 
field coil is partially shielded by the Meissner screening of the sample, reducing the field coil 
- pick-up coil coupling (measured in $o of flux through the pick-up coil per mA of current in 
the field coil). These areas appear dark in the figure. The granular nature of the sample and 
areas of non-superconducting phase are evident. (The lower right area appears most strongly 
superconducting, however this is probably an artifact of topography allowing the SQUID 
closer to the sample.) Fig. 3 also shows an electron backscatter diffraction image, a technique 
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which reveals crystal lattice orientation, of the polished surface of a different piece of the same 
sample. The average grain diameter, in a circle approximation, is 5.1 /xm. 







Fig. 3. Top: susceptibility scan at 4.4 K of the area imaged in Fig. 2. Bottom: electron backscatter 
diffraction scan of a polished surface of a piece of NdFeAsOo.94Fo.o6fi"oii^ the same sample. 

At first glance the mottled background observed in Figs. 2(b)-(e) may resemble the com- 
plex magnetization expected for orbital frustration in a polycrystalline sample. However, to 
the extent visible between vortices, this background is identical in Figs. 2(b)-(e), whereas 
frustration-related moments are polarizeable by cooling in ^T-scale applied fields.^ Instead, 
the background is consistent with an uncancelled in-plane field: as indicated in Fig. 2 we must 
apply ^ —3.3 /xT to cancel the z-axis component of the ambient field, and a comparable 
in-plane component can be expected. This would result in in-plane vortices which would leak 
out near the surface of the inhomogeneous sample. Also, field lines above the sample would be 
deflected upward and downward by the surface inhomogeneity. Both effects would contribute 
to a mottled background signal. 

Polarizable moments would lead to the Wohlleben effect, a bulk paramagnetism against 
the field in which the sample was cooled for fields ^100 /xT, and which has been observed for 
polycrystalhne cuprates.^ To test for the Wohlleben effect in NdFeAsO0.94F0.06, we compare 
the average signal over the sample with the signal beyond the sample edge (in all cases the 
sample was cooled and scanned in the same field). The result, shown in Fig. 4, indicates 
diamagnetism against sub- 100 /iT cooling fields, consistent with an absence of polarizable 
moments. 
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Fig. 4. From each scan, the difference between the mean Bz over the sample and in the corner farthest 
from the sample, against applied field (during cooling and scanning) Ha. 



Qualitative examination of the scans and the absence of the Wohlleben effect indicate an 
absence of orbital frustration in NdFeAsO0.94F0.06- 

In the remainder of the paper, we use established modelling techniques to show quan- 
titatively that the signal that would emerge from orbital frustration would be larger than 
the small mottled background field evident in Fig. 2. The observed vortex widths suggest a 
Josephson penetration depth Xj comparable to the grain size, so spontaneous moments would 
not be well-isolated. Tightly-spaced moments would also tend to align antiferromagnetically, 
further reducing the expected signal at the SQUID. We estimate the expected signal from 
orbital frustration by modelling the grain interfaces as a long 1-D Josephson junction, with 
a single Aj, divided into 0- and vr-junction domains. A 1-D junction is a reasonable approxi- 
mation because the grain size is comparable to the system resolution. In the narrow junction 
limit the phase change across the junction, (/)(x), satisfies a sine-Gordon equation, 

f\2 J. -1 

— = -jj sin((/)(x) + 6>(x)), 

where 9(x) is the position-dependent frustration phase (set here to or tt). 
Two empirical estimates of the typical Xj for this sample are available. 



(1) 



Aj 



(2) 



\ Snedjc 

where d is the magnetic width of the junction and jc its critical current density. (The narrow 
junction limit is d ^ Aj.) d = do + Ai + A2, where do is the actual intergrain spacing and 
Ai and A2 are the penetration depths of the two grains. The grain orientations being random 
these will fall between A^^ and Ac- A^^ has been measured at ^ 200 nm in Tc ^ 50 K Sm- 
and Nd-based samples, ^^~^^ and Xc/Xab ^ 5 has been measured in NdFeAsOo.goFo.io-^^ For 
jc ^ 2000 A/cm^ and d ^ 2 /xm, Aj ^ 4 /xm is obtained. 

The other estimate of Xj comes from the observed vortex widths. The soliton solution, for 
a vortex within the junction, to the sine-Gordon equation is obtained by setting 9(x) to zero 
everywhere. To extend this solution to above the sample we model it as a line of monopole 
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sources an effective height h beneath the pick-up coil, where h is the actual scan height plus 
A, and then integrate Bz over the pick-up coil area, h can be estimated from the susceptibility 
scan shown in Fig. 3: over the sample the field coil - pick-up coil coupling is reduced by ^ 0.16 
$o/niA relative to in vacuum, which would happen with the field and pick-up coils ^ 5 fim 
above a hypothetical A = plane. Setting h — ^ /xm, observed vortex FWHMs of 10-16 jiva 
indicate \j in the range of 1-4 jiva. 

We simulate orbital frustration with a discretized junction 20,000 elements in length, 
divided into domains of mean length L — 40. In each domain 9{x) is set to tt with probability 
P, and zero otherwise. (j){x) is obtained numerically as described in ref. 27. To simulate gradual 
cooling, the system was first solved with \j — 200, then \j was reduced in steps to 10, taking 
the solution from the previous step (with a small perturbation to disrupt unstable solutions) 
as the starting point for the next. L and \j are then scaled to lengths in microns, and the 
solutions are again extended to above the sample by modelling as a line of monopole sources. 

The results of this simulation are shown in Fig. 5, and although it is an approximate 
model the simulation shows that for a wide range of reasonable choices of L, h and \j a signal 
comparable to or larger than the observed background would result: the fixed background is 
very unlikely to be obscuring an orbital frustration signal. The inset compares the expected 
signal distribution for the particular case P — 0.25 and L, /i, \j = 3, 6, 3 /xm, respectively, 
with the observed background (after plane subtraction) in the area indicated in Fig. 2(b). 




Fig. 5. Root-mean-square expected SQUID signal from a frustrated system for various \j\ results of 
the simulation described in the text. Filled circles: P = 0.25; solid line: L, /i = 3, 6 /im, respectively; 
dashed: 3, 4 /im; dotted: 5, 6 /im. Open circles: P = 0.50, L = 3 /im, h = 6 /im. Inset: unshaded 
area: histogram of expected signal (fraction per 1 m<l>o window) for L = 3 /im, h = 6 /im, A j = 3 
/im, P = 0.25; shaded area: observed background in the area indicated in Fig. 2(a) (fraction per 
window divided by 3). 

We have demonstrated that there very likely are no tv phase shifts between tunnelling in 
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different directions in NdFeAsO0.94F0.065 making p and d orders unlikely, s order, vr-shifted or 
not, and s^d order where the d component is small are not ruled out by our result. 
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